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Abstract 

The scalar modes of the geometry induced by dimensional decoupling are investigated. 
In the context of the low energy string effective action, solutions can be found where the 
spatial part of the background geometry is the direct product of two maximally sym- 
metric Euclidean manifolds whose related scale factors evolve at a dual rate so that the 
expanding dimensions first accelerate and then decelerate while the internal dimensions 
always contract. After introducing the perturbative treatment of the inhomogeneities, 
a class of five-dimensional geometries is discussed in detail. Quasi-normal modes of 
the system are derived and the numerical solution for the evolution of the metric in- 
homogeneities shows that the fiuctuations of the internal dimensions provide a term 
that can be interpreted, in analogy with the well-known four-dimensional situation, as 
a non-adiabatic pressure density variation. Implications of this result are discussed with 
particular attention to string cosmological scenarios. 
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1 Introduction and motivation (s) 



It is not unreasonable that internal dimensions of Planckian size can play a role in the early 
stages of the life of the Universe [1] . The question is then how internal dimensions became 
so small in comparison with the external (expanding) dimensions. In the last twenty years 
various sorts of cosmological compactifications have been proposed [2, 3, 4, 5] (see also [1] for 
a review). These mechanisms invoke, in some way or in another, a long phase during which 
the scale factor of the internal manifold shrinks until it stabilizes to the wanted Planckian 
value. There are at least three independent problems associated with a mechanism of 
cosmological conpactification: 

• the nature of the sources driving the contraction of the internal dimensions; 

• the nature of the mechanism stabilizing the radius of the internal manifold; 

• the properties of the inhomogeneities induced by the process of compactification. 

Even higher-dimensional Einstein equations allow for vacuum solutions of Kasner type. 
These solutions may well describe the contraction of the internal dimensions. The caveat 
with such a simple realization of cosmological compactification is the possible occurrence of 
physical singularities both in the external and internal manifold, as repeatedly noticed in 
connection with string cosmological models [5]. 

The properties of metric perturbations induced by a phase of compactification, have 
been partially discussed by various authors. In particular, the Bardeen formalism [7] has 
been generalized to the case of anisotropic manifolds [8, 9] (see also [10]). The evolution 
equations for the Bardeen potentials have been solved in a number of explicit examples 
[11, 12]. Canonical normal modes of the fluctuations have been obtained in specific classes 
of anisotropic manifolds [9, 13]. 

The aim of the present paper is to argue that the contraction of the internal manifold is 
not a sufficient condition for the perturbative treatment of the curvature inhomogeneities 
induced by the process of compactification. In the purely four-dimensional case, curvature 
fluctuations are conserved over wavelengths larger than the Hubble radius in a number of 
relevant cases like, for instance, the one where a single scalar field sources the dynamics 
of the background geometry. If internal dimensions are present, their fluctuations will 
necessarily induce a term analog to the non-adiabatic pressure density variation. Therefore, 
curvature fluctuations will no longer be conserved. This occurrence would not be, by itself, 
problematic. However, a problem may arise if the growth of curvature fluctuations is too 
sharp when the relevant wavelengths are larger than the Hubble radius. 
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To discuss this problem we will use, as a toy model, the low-energy string effective action 
typical of pre-big bang models supplemented by a (non-local) dilaton potential. Within this 
set-up is possible to obtain rather reasonable solutions describing the compactification of 
the internal dimensions and the simultaneous expansion of the external dimensions [14, 15, 
16, 17]. An interesting feature of these solutions is that the external dimensions expands 
in an accelerated fashion and then decelerate (but always expand). The transition from 
acceleration to deceleration is regular in a technical sense. Finally, these solutions can 
be usefully expressed in analytical terms and this aspect helps an unambiguous analytical 
treatment of the evolution of the fluctuations of the geometry. 

The present paper is then organized as follows. In Section 2 the background evolution 
will be specifically described in the case of multidimensional cosmological models charac- 
terized by d expanding and n contracting (internal) dimensions. In Section 3, the salient 
features of the evolution of the fluctuations of the geometry will be derived. Owing to the 
large number of scalar degrees of freedom present in higher-dimensional (anisotropic) man- 
ifolds, the analysis will be then preformed in a five-dimensional geometry. In Section 4 the 
evolution equation of the generalized curvature fluctuations will be presented. Moreover, 
the evolution equations of the system will be recast in a useful first-order form. In Section 
5 the issue of quantum-mechanical normalization of the fluctuations will be discussed and 
quasi-normal modes of the system introduced. Section 6 collects some numerical discussion 
while the concluding remarks will be presented in Section 7. The Appendix contains a series 
of useful technical details on the derivation of the perturbed five-dimensional equations. 



2 Dimensional decoupling 

Multidimensional cosmological models whose line element can be expressed as 

ds^ = GABdx^dx^ = dt^ - a^{t)5ijdx'dx^ - b^it)5abdy^dy\ (2.1) 

will be considered in the present Section. In Eq. (2.1) Gab is the full metric; the indices 
A, B run over the full D = d + n + 1-dimensional space-time; the d-expanding dimensions 
are typically three and the n-contracting dimensions arc generic but the case n = 1 will be 
of particular interest; for notational convenience, in explicit formulae, the indices i,j,k,l... 
run over the d expanding dimensions and the indices a,b,c,d... run over the n internal 
dimensions. 

The background equations stemming from the low-energy string effective action can 
then be written, in natural string units, as [5] 

^ -dH^ -nF^ -V = 0, (2.2) 
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H 

F -- 



Hip, 

Fip, 



lip — ip 



o dV 
nF^ + V - — = Q, 
dip 



(2.3) 
(2.4) 
(2.5) 



where H = a/a and F = h/h; Tp = ip — log \/ —G is the so-caUed shifted dilaton; the overdot 
denotes a derivation with respect to the cosmic time coordinate t. The potential V is only 
function of ^, i.e. V = V{Tp). Eqs. (2.2)-(2.5) can be derived from the following (generally 
covariant) action ^ 



2Af-2 / 



iD-l 



X 



\G\e- 



R + G^'f^Va^'^p^ + V{<p) 



(2.6) 



where 



V = Vie"^ 



1 



\G{w)\e-'^^^^^G'^PdMw)dMw)5U{x) - ip{w)] . (2.7) 



The variation of the action (2.6) with respect to G^,/ and tp leads, after linear com- 
bination, to Eqs. (2.2)-(2.5) valid in the case of homogeneous dilaton field. The general 
(I?-dimensional) derivation of Eqs. (2.2)-(2.5) from Eq. (2.6) has been performed in Ref. 
[18] (see also [19]) and will not be repeated here. The combination of Eqs. (2.2) and (2.5) 
and of Eqs. (2.3)-(2.4) leads, respectively, to the following pair of equations 

^2 ^ IdV 
^ -^+29^ 
{H -F) = {H - F)Tp. 



y = o, 



(2.8) 
(2.9) 



Using Eqs. (2.8) and (2.9) it is rather easy to find solutions of the whole system for 
various sets of potentials. Consider, in particular, the case of exponential potential ViTp) = 
—VQe^f. In this case the solution of the system can be expressed as 

a(t) = {t + \Jt'^ + tl)'^, (2.10) 
h{t) = {t + ^Jf' + tl)-vm, (2.11) 
^ = ip^-^-\u{e + tl). (2.12) 

To satisfy the system (2.2)-(2.5), the free parameters of the solution are subjected to the 
consistency relation e^^^^ ^ ~ Concerning the above solution few comments are in 
order: 

^The choice of natural string units implies that 2\^~^ = 1. 
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• the solutions are defined for t between — oo and +00; 

• the solutions are regular (in a technical sense) since all the curvature invariants are 
regular for every t; 

• the class of solutions expressed by Eqs. (2.10)-(2.12) describe a situation of dimen- 
sional decoupling. In fact while a{t) is always expanding (i. e. d > 0), b{t) is always 
contracting (i.e. 6 < 0) at a dual rate (i.e. b{t) ^ l/a{t)). 

In connection with the last point, it can be easily shown that a{t) expands in an ac- 
celerated way for t < (i.e. a > 0), while the expansion becomes decelerated for t > 
(i.e. d < 0). In fact, taking the limit of Eqs. (2.10) and (2.11) for t — ±00, the leading 
behaviour of the scale factors will be 

1 1 
a{t) {-t) 6(t)~(-t)v^, t<Q; 

a{t)^t'^^ b{t)^t~^^, t>Q. (2.13) 

On the basis of Eqs. (2.13), the case of a five-dimcnsional geometry is particularly 
interesting since, for t +00 the scale factor behaves as a{t) ^ i.e. as if the Universe 
would be dominated by radiation. 

3 Evolution equations of the fluctuations and gauge choices 

In higher-dimensional manifolds characterized by D space-time dimensions, the perturbed 
metric 5Gab has, overall, D{D + l)/2 degrees of freedom. Since D gauge conditions can be 
imposed to fix the origin of the coordinate system, the total number of degrees of freedom 
becomes, after gauge-fixing, D{D — l)/2. If = 4, it is rather obvious to classify fiuctuations 
depending on their transformation properties under three-dimensional rotations. This is 
the choice usually employed in the context of the Bardeen formalism. If D > 4 one can, 
in principle, classify the metric fluctuations according to their transformation properties 
under some larger group of symmetries. If the spatial part of the background geometry 
is decomposed as in Eq. (2.1), it is practical to classify fluctuations with respect to the 
rotations in the expanding dimensions. This was at least a possibility studied in [9, 13] 
(see also [8]). In this approach, as soon as the number of internal dimensions increases 
the number of scalar modes also increases. For practical reasons, the analysis will then be 
conducted in the case of a five-dimensional geometry. Extensions of the formalism to I? > 5 
have been also proposed [9, 13]. 

To derive the evolution of the fluctuations of the geometry it is practical to work directly 
in the conformal time parametrization, i.e. a{r])dr] = dt. Within this parametrization the 
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fluctuations of a five-dimensional geometry can be written as 



0(^00 


= za (p, 






= -a^diP 


- a^Qi, 


SGoy 


= -abC, 




SGij 


= 2a'^tpSij ■ 


- 2a^didjE + a^hij + a^{diWj + djWi), 


SGiy 


= —abdiD 


— obHi, 


SGyy 


= 2b% 





(3.1) 

where the tensor (i.e. hij) and vector (i.e. Wi, Qi and Hi) modes satisfy: 

djhi =0, h^ = 0, 

diW' = 0, diQ' = 0, diH' = 0. (3.2) 

So, hij carries 2 degrees of freedom; Wi, Qi and Hi carry overall 6 degrees of freedom; 
and we also have seven scalars. This implies, as expected, that the parametrization (3.1) 
respects the total number of degrees of freedom of the perturbed five- dimensional metric. 
For infinitesimal coordinate transformations 

x'^^x^ = x^ + e^, (3.3) 

with gauge functions 

= (e°, e\ e^), = (a^eo, -a^Cj, -b^Sy), (3.4) 

the tensor modes of the geometry are invariant while the vectors and scalars transform, 
respectively, as ^ 

Qi — Qi Cii 

Wi = Wi + Ci, 

Hi = Hi- ^dyCi, (3.5) 



and 



(!> = 4> — Cq — Heo 

P = P-e' + eo, 
b I a 
-4 + 7: 



C = C--e'y + -dyeo, 



^The prime denotes a derivation with respect to the conformal time coordinate while the dot denotes, as 
previously reminded, the derivation with respect to the cosmic time coordinate. The relation between the 
two derivatives is simply d/drj = ad/dt. 
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■(/) = '(/> + Heo, 

E = E-e, 
~ h d 

i = ^ + dySy + J^eo, (3.6) 

where the following notation has been used 

ei = die + Ci (3.7) 

such that = 0. On top of the fluctuations of the geometry, the fluctuations of the sources 
should also be taken into account: the gauge transformation for the dilaton fluctuation is 
then 

X = X - 'P'^o (3.8) 

where Sip = x- For instance, the longitudinal gauge [7] can be generalized to the higher- 
dimensional case. For the applications under study it is practical to a adopt the generalized 
uniform dilaton gauge stipulating that 

Wi = 0, D = 0, P = 0, X = 0. (3.9) 

that has been already exploited in various four-dimensional calculations [18, 19]. Notice that 
Eq. (3.9) corresponds to five indipendent conditions since Wi is diver genceless. Different 
gauge choices can be adopted (see, for instance, [9]). However, it has been previously 
shown [19] that these gauge choices are less suitable, already in four dimensions, since they 
lead to divergences in the evolution equations of the fluctuations. The expectation that 
the (generalized) uniform dilaton gauge leads to non-singular evolution equations will be 
confirmed by the present analysis. 

In the uniform dilaton gauge the perturbed equations can be written in explicit terms 
(the details of the derivation are reported in the Appendix). The perturbed form of the 
generally covariant equations can be written as ^ 



sgx + SG 



BM 



G^'^'drZAdNV 



1 
2 



Following common practice, Gmn and denote the background values of the metric and of 

the Christoffel connections, while SGmn and ST'^g are the first-order fluctuations of the corresponding 
quantities. 
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where SQ^ denotes the fluctuation of the Einstein tensor and where 



Equation (3.10) holds in the case of the scalar modes of the geometry since terms like 
57g oc Q'\ relevant only in the analysis of the vector modes [20], have been omitted. Further 
details on the derivation of the evolution equations of the first-order fluctuations are reported 
in the Appendix (see, in particular, Eqs. (A. 10) and (A. 11)). 

The components of Eq. (3.10), computed explicitly in Eqs. (A. 11), give the set of per- 
turbed equations for the evolution of the fluctuations. More specifically, from the perturbed 
(00) and (Oi) components of Eq. (3.10), the Hamiltonian and momentum constraints are, 
respectively, 



and 



+^(^ + J^)dyC + 3^V^V + V^fC + 2V - (^ + H)E' - ^VlE] = 0, (3.12) 

-di[2i;' + r + (^ - m - + n)<f> + ^dyC] = 0. (3.13) 

In Eqs. (3.12)^(3.13) and in the following equations denote the Laplacian with respect 
to the internal coordinate y. 

From the (i = j) component of Eq. (3.10) we have 

^[2jp" + 2H' - 57^2 - - AW - y 
-2^" - f + (^ + 'H)(f>' + 2{Tp' + + {Tp' + 2H- J^)C' 
+V^[E" - (^ + H)E' + V; -f e - </> - ^^lE] 



,2 



+^Vl{2i: -ct>)- ^dy[C' - (^ + H)C] = 0. (3.14) 
Finally the {i / j), (yy), (Oy) and (iy) components of Eq. (3.10) are, respectively, 

did^lE" - (^ + n)E' + ^ + e - - ^VlE] = 0, (3.15) 
-3< + (^+J^)(0' + 3V'') 

+4>[2-^' - ^ V 2T' - 2{T + - '^'HT - m'^ -T"^ - y-^] 

+V2[2V'-0 + ^"-(^' + J^)£;'] = 0, (3.16) 
^V^C - ^a^|3V'' - (^' + J^)0 + 3(7^ - - V\E' + {H- J^)S]| = 0,(3.17) 

2 

^di[c' -{Tp' + 2n- j')c] - '^didy[2^ - 0] = o. (3.18) 
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Concerning the above system of equations the following remarks are in order: 



• if the gradients with respect to y are neglected, the evolution of C decouples; 

• not all the equations are independent; 

• the equations for the zero-modes with respect to y have a slightly more tractable form. 

Let us now prove, in detail, each of the above statements. Since Eqs. (3.17) and (3.18) 
imply, respectively, V^C = and C = {Jp' + 2H — J^)C, indeed the evolution of C decouples. 

Concerning the second remark, from Eq. (3.13) the variable <p can be expressed as 

2i.' + e + iJ^-n)^ a dyC 

(ip' + n) ^ 2h [Tp' + HY ^^-^^^ 

Furthermore, from Eq. (3.15), the variable E" can be expressed as 

E _(c^ +H)S +_V,^; + + _____^_^, (3.20) 

where the dependence upon (p has been eliminated through Eq. (3.19). 

Equations (3.19) and (3.20) can be used to eliminate and E" in Eq. (3.14). The 
remaining dependence on C can be simplified by means of Eq. (3.18). Then, repeated use 
of the background equations imply that Eq. (3.14) is identically satisfied. 

It is here appropriate to recall a useful form of the five-dimensional background equations 
in the conformal time parametrization, namely 

^'=^(^ + 7^) + /(F), (3.21) 

n' - = {n-T){^ + H), (3.22) 

- W? -T'^ = Va^, (3.23) 

where 

depends on the specific form of the potential. Eqs. (3.21)-(3.23) can be easily obtained from 
Eqs. (2.2)-(2.5) and by taking a combination of Eqs. (2.2) and (2.5). Notice, finally, that 
in the case of the potential leading to the class of solutions described in Eqs. (2.10)-(2.12), 
f{V) = Va\ 

Using Eqs. (3.21)-(3.23) and neglecting the gradients with respect to the internal coor- 
dinate, the system of the perturbed equations can be reduced to a more tractable form. In 
particular, Eqs. (3.12) and (3.13) lead, respectively, to 

V2(^ + 2i)) - (^ + n)'^'^E' - Va^ct) + {lp' + J^)(! + 3(^ + U)'^' = 0, (3.25) 
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and 

2V'' + + n)i - + n)<p = 0, (3.26) 

while Eq. (3.14) implies 

-2V'" - + (^' + n)(t>' + 2{Jp' + H)^' + {Tp' + 2n- J^)C' + /(F)0 = 0. (3.27) 

Finally, Eqs. (3.15) and (3.16) become, respectively, 

E" - {Tp' + n)E' + V' + e- </> = 0, (3.28) 

-3^ + (^' + J^)((/'' + 3V'') - f{V)cP + V2(V' - - (-^ - W)V2£' = 0. (3.29) 

4 Generalized curvature fluctuations 

For the purposes of the subsequent investigation, it is practical to define the variable 

A = V + |. (4.1) 



In terms of A, Eq. (3.19) becomes 



(^' + H) ip' + H 



Inserting Eq. (4.2) into Eq. (3.25) and using Eqs. (3.21)-(3.23) the generalized evolution 
equation for the curvature parturbations can be obtained: 

^ (^^+3w-2^)(^^+w) _ (.F-?^)(3?^^ + ^^-^^^) ^ 

+ V^A+ „<^ + :<'^^ VV. (4.3) 



(^ + 37^)2 + 2j^2 (^/ + 37^)2 + 2J^2 

Clearly, in the limit ^ and ^ ^ (i.e. X ip), Eq. (4.3) reduces to 

which is the typical form for the evolution of curvature perturbations in the uniform dilaton 
gauge [19]. Eq. (4.4) suggests that if ip' = {Tp' + 3?^) / (as in the class of models 
discussed in [19]), then V' — for long wavelengths larger than the Hubble radius. If 
internal dimensions are dynamical, Eq. (4.3) generalizes the evolution of the curvature 
perturbations. 
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Using the definition (4.1) and inserting Eq. (4.2) into Eq. (3.28), Eq. (4.3) allows to 
eliminate A', and the resulting equation will be 



E" 



(^ + H)E' + 



2{Tp' + n) 



W + 37^)2 + 2J^] 



A + 



+ 



(4.5) 



iip' + mf + 2j^ ^ 2[{Tp' + my + 2j^] 

Since the background evolution has a simple analytical form in t (but not in ??), Eqs. (4.3) 
and (4.5) can be written, for future simplicity, in the cosmic time coordinate and in terms 
of the auxiliary variables 

£ = c?E, V = t (4.6) 

whose convenience will become clear in the context of the numerical analysis. The result of 
this simple algebra will be, in Fourier space, 



{Tp + H){Tp + m - 2F) 
2[(^ + 3i7)2 + 2F2] 



{Tp + my + 2F2 

2{Tp + H) ,.2, 



(^ + 3iJ)2 + 2F2" ' (^ + 3ii")2 + 2F2" 
where the physical frequency of the fluctuation, i.e. u = k/a has been introduced. With 
similar manipulations, the evolution equation for £ will be instead: 



Tp + 2H 



(^ + 3i?)2 + 2F2, 



311 - 2F 



{Tp + 3i?)2 + 2F2 



[^^ + H{lhH + 12^) + 2^2 _ 6ir(^ + H)] 
2[(^ + 3iJ)2 + 2F2] 



ik- 



1 - 



4a;2 



[Tp + mf + 2F'^, 



Afc. (4.8) 



Equations (4.7) and (4.8) define a first-order linear differential system (with time- 
dependent coefficients) in the variables, £k, Xk, Vk and ^k- Having already specified the 
equation for ^j., i.e. 

ik = Vk, (4.9) 

we just need a further relation determining Vk- After some algebra (which will be swiftly 
mentioned in the following section in connection with the problem of the canonical normal 
modes) the wanted evolution equation turns out to be 

F{Tp + 3H - 2F) 



-{ 



Vk = <^ + 2V-^ 



{Tp + 3H + F) [{Tp + 3Hf + 2F'^] 
{Tp + 3H){Tp + 3H - 2F) 



Vk 



{Tp + 3H + FY[{Tp + 3HY + 2F2] 



+ 



^^^ {F-H){Tp + 3H^- 2F) {3H^ + ^2 + Hp) ^ 



{Tp + 3H + Ff [{Tp + 3Hf + 2F2] 



ik 



-4V 



{Tp + 3H- 2F) [F2 + 3ff 2 + Hp] 
{Tp + 3H + Ff [{Tp + 3Hf + 2F2] 



2 c 
LO tk, 



(4.10) 
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where the expHcit form of the potential mentioned before Eq. (2.10) has been used. 

Concerning Eqs. (4.7), (4.8) and (4.10) it is appropriate to recall that the time- 
dependent coefficients are always regular for any finite value of the cosmic time coordinate 
as it can be argued by noticing that, in the denominators, the only combinations that appear 
are 

(i^ + SH + F), {^ + 3Hf + 2F^. (4.11) 

It is clear that the second expression in Eq. (4.11) never goes to zero since it is the sum 
of two positive definite quantities. Moreover, the first quantity, corresponding to ip is also 
positive definite for the class of backgrounds obtained by setting d = 3 and n = 1 in Eqs. 
(2.10) and (2.11). 

Equations (4.7), (4.8) and (4.10) can then be numerically integrated once the initial 
conditions for the fluctuations are set in the limit t ^ —oo. Consequently, it is mandatory 
to determine which are the correct canonical normal modes that have to be used to enforce, 
for instance, quantum mechanical initial conditions for the fluctuations. 

5 Quasi-normal modes 

The fluctuations of the geometry have to be normalized in the limit 77 — —00 (or t —00) 
and, for this purpose, the canonical normal modes have to be determined. In the system 
under consideration, when t — ^ —00 the potential term becomes negligible with respect to 
tp^, i.e. y ^ (^^. In this situation, two normal modes of the system can be identified. They 
will be called quasi-normal modes since, for t ^ to (when the potential dominates) the two 
modes mix non-trivially. 

As far as the evolution of the system is concerned, it is necessary to discuss directly the 
system derived in Eqs. (4.7), (4.8) and (4.10). The reason why this choice is not arbitrary 
stems from the particular form of the equations characterizing the normal modes of the 
system. These equations contain time-dependent coefficients whose denominators include 
powers of either {p + H) or (Tp + F). Now, these terms lead, generically to poles in the 
time-dependent coefficients of the system. For instance, in the specific case of the class of 
solutions (2.10)-(2.12) the zeros of {jp + H) and (jp + F) are, respectively, in {d + n— 1)"*^/^ 
and —{d + n — 1)^^/^. For the specific case of the numerical analysis to be presented in the 
following section this would imply poles (of various degree) for either 1/^/3 of -l/\/3. 

To prove the statements of the previous two paragraphs, let us derive the evolution 
equations for the normal modes of the system first by looking at the evolution equations 
and then by looking directly at the action perturbed to second order in the amplitude of 
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the fluctuations. Summing up Eqs. (3.25) and (3.27) and eliminating (p by means of Eq. 
(3.28) the following equation can be derived: 

V^X-X" + 3ilp' + H)X' + 9{V)(j) = - '^A^ [^'" - + n)E"- {ip" +n')E'- V^E'] , (5.1) 



9iV) = ^{^-W] (5.2) 



where Eq. (4.1) has been used and where 

4 

is a new function depending on the specific form of the potential. Notice that g{V) vanishes 
for the class of solutions described in Section 2. Taking the difference of Eq. (3.29) and 
(3.27), the resulting equation is, in some way, similar to Eq. (5.1): 

ie - O - v2(c - V) + i^'ijp' + ^-m)- e'(^ + m- 2t) = 

-(J^ - n)[E"' -{Tp' + H)E" - {Tp" + n')E' - V'^E']. (5.3) 

Hence, combining Eqs. (5.1) and (5.3) to get rid of the terms containing the derivatives of 
E, the following simple relation appears 

W^-+(^+7^)e-+^(^)g(F)0 = ^( ^^+_f^^'^ )[v2A-A-+(^+7^)y]. (5.4) 

This equation can be used either to eliminate the Laplacian of ^ in favour of a Laplacian of 
A or vice-versa. 

To get a quasi-decoupled equation for A, the first step is to take the conformal time 
derivative of Eq. (4.3). Terms like V'^E" will naturally appear and their presence can 
be eliminated by means of Eq. (3.20) or (4.5). As a result of this manipulation terms 
proportional to V^^ and to V'^E' will arise. They can both be eliminated by the use of Eqs. 
(5.4) and (4.3). The final result for the evolution equation of A reads then 



A" 



1 + -a 



3 odV 1 



2 d(p{(p' + 3n + 



V^-A' 



(^ + W) + 



(^' + 37^ + J^)2 ^^^'^^ 2((^' + 37^ + J^)2'* 

+^ (jp' + nw + sn + j-r ^ ^ ^ 

As previously remarked, in the limit rj — ^ — oo the potential terms are negligible with 
respect to the other terms of comparable dimensions such as Tp''^, Ti?, and their mutual 
combinations. Consequently, the evolution equation for A becomes, in this limit 

A" + 2^A' - V^A = 0, (5.6) 
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where _ 

.,^-,Jl^t4^. (5.7, 
V 3 va if +n 

Notice, as anticipated, that if we ought to use Eq. (5.5) for expHcit numerical integration 
over the whole range of 7], poles arise for {Tp' + 7i). This confirms that second-order equa- 
tions such as (5.6) are mandatory for the determination of the initial conditions but not 
appropriate for the description of the evolution of the fluctuations. 

The same type of algebra leading to Eq. (5.6) can be applied to find the second (asymp- 
totic) normal mode of the system. Indeed, the variable 

also obeys the following equation 

e" -{Tp' + n)9' - v'^e = 

^^"^ [2(^ - n)\' - {2H + 7^)9'] 



Va'^{2n + F) 



[^^ -n^- 2Va^\e, (5.9) 



{Tp' + m + T){Tp' + nY 

whose specific form becomes, in the limit rj — ^ — oo, 

z' 

Q" + 2^9' - V^9 = 0. (5.10) 

Zg 

Hence the canonical normal modes are 

V = ZgX, w = Zs9 (5-11) 
and their Lagrangian density can be written as 



(5.12) 



It is relevant to stress that Eq. (5.12) holds, strictly speaking, only in the limit of vanish- 
ing potential (i.e. rj — ^ — oo) and agrees with the Lagrangian for the normal modes of a 
dimensionally reduced geometry discussed in [13]. As a side remark, we ought to conclude 
the present section by noticing that, using Eq. (5.5) into Eq. (5.4), the equation 



{Tp' + 3n + J^) [{Tp' + mf + 2J^^] 

, {Tp' + mw + m-2T) , 

{Tp' + m + j')[{Tp' + mf + 2j^Y 
2 {^ - nw + m- 2T){m^ + :f^ + W) ^ _ ^ 

(^' + 37^ + J^)2[(^' + 37^)2 + 2jr2] ? ^ ? 

{^' + 3n + m{^' + m' + 2:FY ^ ^ 

can be obtained. This is the equation already anticipated in Section 4 (see Eq. (4.10)). 
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6 Numerical analysis 



To solve numerically the evolution it is essential to have a set of evolution equations whose 
time-dependent coefficients are all non-singular. The relevant set of equations follows from 
Eqs. (4.7), (4.8) and (4.10) recalling the remark related to Eq. (4.11). 

To solve this system numerically, quantum mechanical initial conditions have to be 
imposed for t — oo. According to the results of the previous section, the Fourier modes 
of 6 and A will be normalized as 



1 



-i J w{t')dt' 

2uJ{^/aZs) 

1 -if uj{t')dt' 



(6.1) 



Prom Eq. (6.1), the asymptotic expressions of Ok and of can be deduced by taking the 
derivative of the two expressions with respect to t. Having determined the values of 0^, A^ 
and of their derivatives it is possible to set the initial conditions for the other variables. The 
function Afc directly appears in the system. Recalling the explicit expressions of Afc and Ok 
in terms of and ■i/'fe, i-e. Eqs. (4.1) and (5.8), the initial conditions for and Vk can be 
easily obtained: 



ik 



Vk = l 



Tp + m + F 

Tp + H 

ip + m + F 



Gk 



ip + m -2F 

Tp + H 
Tp + m ~2F 



+ 



Tp + H J V Tp + H 
2 V[2{F-H)Xk-{2H + F)0k 
3 i^ + Hf 



As 



(6.2) 



The last quantity to be dtermined is the initial condition for Sk whose expression, in terms 
of the initial conditions of the other variables, is: 



■ 



i<^ + 3H 



2uj'^{Tp + H) 
3Hf + 2F^ 



2F) _{Fj~H)^H^ + F^-fl 



uj^iTp + ny 



LU^{Tp + H)^ 



(6.3) 



The initial time of integration, ti, is fixed in such a way that the given physical frequency 
is larger than the Hubble rate, i.e. the corresponding wavelength is smaller than the Hubble 
radius. In the case of the class of solutions introduced in Eqs. (2.10)-(2.12) this requirement 
implies, for a given Fourier mode k, ti < — (2av)~^/^ to where k = Hq. Thus, as k becomes 
smaller, the initial integration time becomes larger in absolute value. 
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Figure 1: The evolution of the fluctuations is ihustrated in the case of a five-dimensional 
background belonging to the class (2.10)-(2.12). We took ipQ ~ —40. As previously intro- 
duced K = ktQ. 

Consider now the five-dimensional background geometry obtained from Eqs. (2.10)~ 
(2.12) by setting d = 3 and n = 1. In this case Eqs. (4.7), (4.8) and (4.10) can be 
integrated by setting quantum mechanical initial conditions as discussed in the previous 
paragraphs of the present section. The result of the numerical integration is illustrated in 
Fig. 1 where the logarithm (in ten basis) of spectrum of A, i.e. 



is reported. In analogy with Eq. (6.4) we will also define, in similar terms, the spectrum of 
^, i.e. 6^ and so on. 

It is clear from Fig. 1 as well as from the equivalent plots obtained for 6^ and Ss that, 
in spite of the initial conditions, for t > the following approximate relation holds 



where a is a constant and it is numerically determined to be 4.9; the symbol of similarity in 
Eq. (6.5) means that the relation holds up to constant factors of order 1. This occurrence 
signals that there is a kind of late-time attractor in Eqs. (4.7), (4.8) and (4.10). According 
to Eqs. (6.5), the plots for the other power spectra will be quantitatively and qualitatively 
similar to the one of A. The sharp (exponential) amplification of the fiuctuations, is rather 
dangerous for the validity of the perturbative expansion, as it can be argued from Fig. 1. 



(6.4) 



(6.5) 
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Figure 2: The evolution of curvature fluctuations is illustrated in the case where F ^ for 
l>t\ with l\ = 10. 

To gain analytical understanding of the behaviour reported in Fig. 1 let us write the 
system provided by Eqs. (4.7), (4.8) and (4.10) in the limit of large (positive) times, i.e. 
t ^ to. Denoting with r = t/to and with k = kto the system becomes: 

■ /2 , 17k2\ / 1\ /3 10 \^ /8 ^ 2k2 \ 

£k ~ -K 7T <Sfc + 2r H ]Vk+ { ^ kfc + -1^ T H 1 Kfc, 

V3 18r2y " \ 3tJ ^ V2 3t^ \3 9r J 

ik ^ Vk, 

In Eq. (6.6) we kept all the leading terms in the limit r ^ 1 as well as the first subleading 
correction. The system (6.6) can be solved in the limit k — > 0, i.e. when all the modes have 
wavelengths larger than the Hubble radius at tQ. In this limit the system becomes 

4-24^^, Vk^ik, Sk-^Ck, Afc~-^Pfc. (6.7) 
The system (6.7) can be easily solved and the result is that an exponential mode 

Ck{t) ^ Uk) e^^^ Xk{t) -Ve^oik) e^v^^ £k{t) ^ ^ ^k) e'^^\ (6.8) 

may arise. Equation (6.8) is the rationale for the result of Eq. (6.5) and, indeed, 2\/6 ~ 
4.8989. In Fig. 1, the dot-dashed line is nothing but the slope of the logarithm (in ten 
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Figure 3: The same analysis discussed in Fig. 2 is here repeated in the case t\ = to- 

basis) of the function e^^(*/*o)^ This example shows that the fluctuations of the internal 
dimensions provide an effective non-adiabatic pressure variation whose role is to provide a 
source term for curvature perturbations. As it is evident from this example, the continued 
contraction of the internal scale factor is not sufficient for a perturbative description of the 
fluctuations. In other words, it is the fact that F does not vanish for t > Iq that triggers 
the exponential growth of the curvature fluctuations. 

To stress this point, let us now consider the case when the solution is exactly the 
one previously studied but only for t < ti. For t > ti the Universe is described by a 
five-dimensional geometry where, however, b ~ constant and F is driven dynamically to 
zero. Under this assumption let us now solve the evolution equations of the fluctuations, 
i.e. Eqs. (4.7), (4.8) and (4.10). Of course, as discussed above in the present and in 
the previous section, quantum mechanical normalization will be enforced for t — > — cxd. In 
Fig. 2 the result of such an exercise is reported in the case ti = W to. It is clear that the 
dangerous exponential amplification is absent and that the treatment of the inhomogeneities 
is perturbative while the relevant wavelengths are larger than the Hubble radius. If b{t) 
freezes at an earlier epoch, the amplitude of the curvature fluctuations will be smaller. This 
observation is corroborated by Fig. 3 where the same quantity illustrated in Fig. 2 is 
computed in the case ti — to- Since ti is smaller, F ^ earlier than in the case of Fig. 2. 
This implies that the (short) phase of exponential amplification present in the example of 
Fig. 2 is now absent. 
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7 Concluding remarks 



In the present paper the possible occurrence of dimensional decoupling has been scrutinized 
in light of the possible space-time fluctuations of the internal dimensions. Dimensional 
decoupling arises in various classes of higher-dimensional (anisotropic) manifolds character- 
ized by the expansion of the observable dimensions and by the contraction of the internal 
dimensions. In the present paper, dimensional decoupling has been described by means of 
a class of solutions of the low-energy string effective action. The virtues of this description 
are that the background solutions arc always regular. Therefore, a reasonably accurate 
description of the evolution of the inhomogeneities can be developed. 

Various results have been obtained: 

• the full system of the equations for the perturbations has been derived in the string 
frame description; 

• the evolution equations of the zero-modes have been diagonalized in terms of a set of 
variables that are non-singular throughout the evolution of the background geometry; 

• quasi-normal modes of the system have found; 

• a numerical analysis of the appropriate (first-order) evolution equations has been 
presented. 

The obtained results show that the contraction of the internal dimensions is not a suf- 
ficient requirement for the perturbative description of the inhomogeneities. If the scale 
factor of the internal dimensions is not fixed to a constant (asymptotic) value, then the 
evolution of the internal fluctuations triggers the presence of a source term for the evolution 
of curvature fluctuations. This term may be interpreted as the analog of a non-adiabatic 
pressure-density fluctuations whose role, in four space-time dimensions, is known to imply 
a growth of curvature fluctuations. A possible problem arises, in this context, if the con- 
traction of the internal manifold lasts for a long time. In this case the numerical results 
show that the evolution equations of the fluctuations may develop some kind of dynami- 
cal attractor where the curvature fluctuations are exponentially amplified in spite of the 
initial conditions. The numerical analysis can be corroborated, in the specific case under 
consideration, by a precise analytical understanding. This potential problem of dimensional 
decoupling can be avoided provided the internal scale factor sets, sufficiently early, to a con- 
stant value. In this case the evolution of the fluctuations may still be perturbative during 
the whole evolution. 
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The results reported in the present paper suggest that reahstic models of multidimen- 
sional pre-big bang dynamics can be constructed. However, unlike the attitude taken so 
far, care must be taken in dealing with the fluctuations of the internal dimensions that can 
induce dangerous late-time attr actors in the evolution equations of the fluctuations. 
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A Technical details 



In the case of the geometry given in Eq. (2.1), assuming that both the external and internal 
manifolds are spatially flat, the Christoffel connections can be written as 

^00 = ^) ^ij = 'Hdij, Tqj = Hdl, 

r'ib = ^^Sab, ri, = TSl (a.i) 

while the components of the Ricci tensor are: 

Ri = -^[n' + id- l)H^ + nHr\5l 

RI = -^[F' + {d- 1)HJ^ + nJ^^]Sl (A.2) 

From Eq. (A.2), the Ricci scalar is then 

R = -i.[2dn' + 2nJ^' + n(n + l)J^'^ + did - 1)H^ + 2n(d - l)nJ^]. (A.3) 

In the uniform dilaton gauge (see Eq. (3.9)) the non-vanishing components of the perturbed 
metric (3.1) become, 

SGoo = 2a^(l), 

6Goy = —abC, 

SGij = 2a^ip5ij - 2a^didjE, 

SGyy = 2b^^ (A.4) 

and 



5G^y 



2(t> 

c_ 



62- 

The fluctuations of the Christoffel connections become then 



00 



SG^y = (A.5) 



20 



Sr'^j = -[^|J'6ij + 2W((/. + i:)]5ij + didj[E' + 2nE], 

5rjo = -^'51 + did^E' 

}? ■ 



2 



2a 

6Tly = -dy^5i+dydid^E, 

5T'lj = d^il)5ij - djil^S^ - dii;S^ + didjd'^E. (A.6) 

The first order fluctuations of the Ricci tensor with mixed indices are: 

5i?o = lv2[0 - E" - HE'] + lv^0 + 

\dy{C' + TC) + ^[3V" + e" + {m^)' + (2J!^ - U)0\ 
ah a^ 

1 



■^^{m! + 2T' - 2Hr+2r'^), 

= -\d'\2ip' + + ^{T -H) + {2H + r)<i> + ;^C], 

(X ^CL 

54 = ^510^" + [2H' + + 2Hr]<t> + H<t>' 

2 

+{hH + j^)^' + H^' - vV - ^v^V' - T^v^^;' + ^ndyC} 



+did^[-E" - {2n + j=)E' + - ^ - V + ^v^^;], 



2 
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SRl = ^[3n' - 2HJ^] + ^v^c 
ab 2ab 

-^dy{3i;' + 3n(p + 3in - - v\E' + {n- j^)e]}, 

SR'i = ^di[C' + {n + 2J^)C] - ^didy[2i; - </,], (A.7) 
while the fluctuation of the Ricci scalar becomes: 

5R = ^1 V'[</. -^-2^-E"- {Sn + J^)E' + ^VlE] 

+{6n' + 2T' + + + 2j''^)4> 

- 3V') + l^dy[C' + m + ^)]. (A.8) 



Defining now 



Ga = SRl - l6l6R, (A.9) 
the scalar fluctuation of the Einstein tensor leads to the following components: 

Sg^o = ^V^[e + 24^ + {2H + J')E' - ^VlE] 

-\ndyC + + \{-3{2'H + T)^' - - m^n + 

ab b'^ ^ 

6gi = J^5^' |-2V;" - ^" - 2[n^ + + nj^+j^' + 2n']^ 
-(j^ + 2n)<i)' - 2{T + 2n)^' -(n + 2j^)^' 

+V\E" + (27^ + T)E' + V + ^ - <^ - -^V^^;] I 

+|v^(2V' - 0) - l5,[C' + (27^ + ^)]<J^ 

-^d.^^{E" + (2W + + V; + C - - ^V^^], 

= ^j-SV-" - - - 6{H' + H^)(f) 
+V^[2i}-(P + E" + 3nE']y (A. 10) 

where, clearly, SQ^ = SR^ SQ^ = SRl, SQ^ = SR^. 

Finally, the general form of the perturbed equations can be usefully expressed in a 
compact way as 

a^Sg^ + (t>[2J^'fi' - ip'^ + qW] 
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+3</?V' + <P'S! - ^'V'^E' + -^'dyC = 0, 



( ffl BV 

a'^sgy + 4>[2ip" - ip'^ + m<p' - y ^] + 0V' + 3</^V' - "^'v'^;' = o, 
a'sg'o + \c\-^" + + + y ^] + ^'^^'<^ = 0' 

a^^af - ^9^'9,C = 0, 

o^^Q'q + V^'^V = 0. (A.ll) 

Equations (A.IO) together with Eqs. (A.ll) can be used to derive the evolution equations 
of the fluctuations discussed in the bulk of the present paper. 
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